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Abstract 

Ever since their discovery the regular satellites of Jupiter and Saturn have held out the 
promise of providing an independent set of observations with which to test theories of planet 
formation. Yet elucidating their origin has proven elusive. Here we show that lapetus can 
serve to discriminate between satellite formation models. Its accretion history can be under- 
stood in terms of a two-component gaseous subnebula, with a relatively dense inner region, 
and an extended tail out to the location of the irregular satellites, as in the SEMM model 
of Mosqueira and Estrada (2003a,b). Following giant planet formation, planetesimals in the 
feeding zone of Jupiter and Saturn become dynamically excited, and undergo a collisional 
cascade. Ablation and capture of planetesimal fragments crossing the gaseous circumplane- 
tary disks delivers enough collisional rubble to account for the mass budgets of the regular 
satellites of Jupiter and Saturn. This process can result in rock/ice fractionation as long as 
the make up of the population of disk crossers is non-homogeneous, thus offering a natural 
explanation for the marked compositional differences between outer solar nebula objects and 
those that accreted in the subnebulae of the giant planets. For a given size, icy objects 
are easier to capture and to ablate, likely resulting in an overall enrichment of ice in the 
subnebula. Furthermore, capture and ablation of rocky fragments become inefficient far 
from the planet for two reasons: the gas surface density of the subnebula is taken to drop 
outside the centrifugal radius, and the velocity of interlopers decreases with distance from 
the planet. Thus, rocky objects crossing the outer disks of Jupiter and Saturn never reach 
a temperature high enough to ablate either due to melting or vaporization, and capture is 
also greatly diminished there. In contrast, icy objects crossing the outer disks of each planet 
ablate due to the melting and vaporization of water-ice. Consequently, our model leads to an 
enhancement of the ice content of lapetus, and to a lesser degree those of Titan, Callisto and 
Ganymede, and accounts for the (non-stochastic) compositions of these large, low-porosity 
outer regular satellites of Jupiter and Saturn. For this to work, the primordial population 
of planetesimals in the Jupiter-Saturn region must be partially differentiated, so that the 
ensuing colhsional cascade produces an icy population of > 1 m size fragments to be ablated 



during subnebula crossing. We argue this is likely because the first generation of solar neb- 
ula ~ 10 km planetesimals in the Jupiter-Saturn region incorporated significant quantities 
of ^^Al. This is the first study successfully to provide a direct connection between nebula 
planetesimals and subnebulae mixtures with quantifiable and observable consequences for 
the bulk properties of the regular satellites of Jupiter and Saturn, and the only explanation 
presently available for lapetus' low density and ice-rich composition. 
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1 Formation of the Regular Satellites of Giant Planets 



The regular satellites orbit in the prograde direction, meaning the same sense of rotation as 
the spin of the primary, and lie close to the Laplace sm^facc, with the exception of Saturn's 
lapetus (whose orbit is nearly circular with an eccentricity of e = 0.028 but inclined i = 
7.5° with respect to the Laplace surface). Broadly speaking, the interpretation of these 
observations is straightforward: unhke the irregular satellites, which tend to have inclined 
and distant orbits and represent a population of captured objects, the regular satellites 
form in circumplanetary disks of gas and solids. Notably, the regular satellite systems of 
Jupiter and Saturn share a number of similarities, including the mass ratio of the largest 
satellites to the primary, the specific angular momentum, and bulk compositions. Yet, the 
differences are also striking: a trend of decreasing density with radial distance is apparent 
for the Jovian but not the Kronian regular satellites; and the mass ratio between Ganymede 
and the other Galilean satellites is not nearly as extreme as that between Titan and its 
neighbors. Indeed, Titan's isolation, the compositional diversity of Saturn's inner moons 
as well as its ring system may attest to titanic clashes between the well-behaved residents 
proper and unruly foreign hordes - interlopers from the outer solar system wreaking-havoc 
upon the planet's orderly retinue. In its remote outpost, lapetus, a survivor of the onslaught, 
icily bears witness to that early violent epoch, as its battle-scarred surface records. lapetus' 
extreme albedo contrast, large separation from Saturn and other Saturnian moons, its low 
eccentricity yet significant inclination, its frozen-in shape and equatorial ridge, its ancient, 
cratered surface, and icy composition all combine to make it a Rosetta moon. It is to this 
mysterious, two-faced sentinel of the outer solar system that we primarily turn our attention. 

Observational evidence indicates that the largest KBOs are of different composition than 
the regular satellites of the giant planets. Triton (p = 2.061 ± 0.007 g cm~^; Person et al. 
2006), Eris (p = 2.3±0.3 g cm'^ Brown and Schaller 2007), and Pluto-Charon (p = 1.94±0.09 
g cm~^; Buie et al. 2006) have densities that imply a rock/water-ice ratio of approximately 
70/30 by mass. These high densities have long been interpreted to result from accretion 
in the outer solar system given the presence of abundant, oxygen-sequestering uncondensed 
GO (Prinn and Fegley 1981; McKinnon et al. 1997). In contrast, Ganymede (p = 1.94 g 
cm~^) and Callisto (p = 1.83 g cm~^) each consist of ~ 50/50 rock/water-ice by mass, with 
Ganymede being slightly more rock-rich (Sohl et al. 2002). Similarly, Titan (p = 1.88 g 
cm~^) is also half-rock/half-ice by mass. The simplest explanation of these observations is 
that the subnebulae of the giant planet were enriched in water-ice compared to the outer 
solar nebula, although this interpretation is partly clouded by uncertainties in our present 
understanding of solar abundances (e.g., Wong et al. 2008). 

Further insight can be sought from other regular satellites. The inner satellites of Jupiter, 
lo and Europa, may have lost volatiles either due to the temperature gradient in the sub- 
nebula (Pollack et al. 1976; Lunine and Stevenson 1982), coUisional processes involving dif- 
ferentiated objects, and/or the Laplace resonance^. Switching over to Saturn, the observed 
densities of the medium-sized regular satellites are not compatible with solar composition 
(Wong et al. 2008). It has long been argued that these observations are consistent with 

^Amalthea may have incorporated water-ice (Anderson et al. 2005), but it is unclear how much. Amalthea 
may be too small and close to the planet to provide a useful constraint on satellite formation theories. 
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accretion in a reducing subnebula with efficient conversion of CO to CH4 (as opposed to 
inhibited conversion in the solar nebula; Prinn and Fegley 1981). However, sateUite forma- 
tion models do not provide an environment that can support such chemistry (Estrada et al. 
2009). Also, closc-in to Saturn coUisional processes and/or resonances (at least in the case 
of Enceladus) appear to have resulted in a stochastic component^. 

Yet, lapetus would not be affected by these processes, and it is large enough to provide a 
sample of the composition of the outer disk of Saturn. lapetus is different from other regular 
satellites largely by virtue of its distance from the planet: it may have survived the early 
bombardment because impact speeds at its location are significantly slower than those closer 
in to Saturn, where the gravitational well leads to hypervelocity impacts; its frozen-in shape 
and unique ridge are likely associated with its location far from Saturn (Castillo- Rogez et al. 
2007); its significant inclination is likely connected with the deviation of the Laplace plane 
away from the planet's equator plane at its location (Ward 1981 but note that this scenario 
requires a subnebula dispersal time of < 10^ yr; Tremaine et al. 2009); its large albedo 
contrast is likely ultimately due to the deposition of exogenous material on the surface of 
this tidally-locked regular satellite (Burns et al. 2009; Buratti et al. 2009); and, as we show 
here, its composition too is tied to its location in the Kronian subnebula. These properties 
taken together with the difficulty inherent in the capture of a such a large body into a low- 
eccentricity orbit provide a compelling case for in situ accretion in a circumplanetary disk 
of gas and dust^. 

To establish a connection to the outer Kronian subnebula we must first constrain lapetus' 
porosity. Given its mean radius and mean density (r = 736.0 ± 2 km and p = 1.081 ± 18 g 
cm~^; Thomas et al. 2006) the pressure at the core (P = (27r/3)p^Gr^ ~ 100 MPa) is high 
enough that porosity should be small there (though it could be larger in the outer layers of 
the satellite). Experimental results indicate a reduction in porosity in cold, granular water 
ice over a pressure range of ~ 1 — 150 MPa, with porosities of the order of ~ 10% for ~ 100 
MPa pressures (Durham ct al. 2005). Heating, annealing and creep-driven water-ice flow can 
result in pore collapse and decrease the porosity even more. Models in which lapetus de-spins 
due to the presence of short-lived radioactive isotopes (SLRl), leading to the formation of 
a ridge, imply warm-ice and low porosities (Castillo-Rogez et al. 2007; Melosh and Nimmo 
2009). Although the rheology of icy satellites is poorly understood, the interpretation of 
lapetus' fossil shape as "frozen-in" due to lithospheric hardening as the satellite was in the 
process of despinning is persuasive, and at odds with with a high-porosity object. Note that 
weak radiogenic heating due to long-lived radioactive isotopes (LLRl) alone would be enough 
to collapse the pores except in a surface layer tens of kilometers thick (Castillo-Rogez et al. 
2007). An lapetus' model with a porosity of ~ 0.3 in a clcan-icc upper ~ 50 km (roughly 
the surface layer thickness of the Castillo-Rogez et al. (2007) LLRl model; cf their Fig. 6) 
of the satellite would result in a rock ratio of ~ 30% by mass, which is still rock-depleted 
compared to Ganymede, Titan and Callisto. Furthermore, even if it were possible to despin 
a cold lapetus, the presence of the ridge is hkely to indicate ice mobility early on. Thus, we 

^See Mosqueira and Estrada 2003b for an alternative explanation for the ice-rich composition of close-in, 
mid-size Saturnian regular satellites. 

^Elsewhere (Mosqueira and Estrada 2005), we discuss an intriguing coUisional-scattering scenario for the 
origin of lapetus; however, we favor the present model. 
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argue that lapetus is an icy satellite with rock/water-ice fraction of ~ 20% by mass, making 
its composition incompatible with solar mixtures (Wong et al. 2008). 

We stress that lapetus' low density relative to Titan, Ganymede and Callisto cannot 
be explained by a snowline argument since all four satellites are taken to accrete outside 
the snowline (Pollack 1976; Lunine and Stevenson 1982). The issue here is the depletion 
of rock, as opposed to the presence of volatiles, in lapetus. In this regard, it is instructive 
to compare lapetus to irregular Phoebe. The Cassini flyby of Phoebe yielded a density of 
1.630 ± 0.046 g cm~^ (Jacobson et al. 2006). If we allow for a moderate porosity, which 
is reasonable given its 106.60 ± 1.00 km size, this density corresponds to a rock/water- ice 
fraction similar to that of large KBOs. The contrast between lapetus and Phoebe reinforces 
the interpretation of Phoebe as a captured moon from the outer solar nebula (Johnson and 
Lunine 2005). Indeed, the rock/water-ice fraction for Phoebe may be larger than those of 
Ganymede, Callisto and Titan. Allowing for significant porosity in the case of Phoebe would 
accentuate the compositional contrast with lapetus and other regular satellites. 

In the nucleated instability model of planet formation (e.g., Bodenheimer and Pollack 
1986; Pollack et al. 1996) a core must first form by accretion of planctcsimals. In this mode of 
planet formation most of the mass of solids resides in planetesimals of size ~ 10 km (Wetherill 
and Stewart 1993). Planetesimals are also needed to explain the observations of the Oort 
cloud and the scattered belt (e.g., Charnoz and MorbideUi 2007), to power planet migration 
in the Nice model (e.g., MorbideUi et al. 2008), and to explain the volatile enhancement in 
giant planet atmospheres by planctesimal trapping and delivery (e.g., Owen and Encrenaz 
2003). Here we seek to provide a direct physical link between planetesimals in the solar nebula 
and the circumplanetary disks of giant planets, and to account for the source of the solids 
that ultimately led to the formation of the regular satellites of the giant planets. Our aim is 
to study mass delivery by ablation of planetesimal fragments crossing the circumplanetary 
disk (to be followed by re-condensation and satellite accretion). In particular, ablation can 
result in fractionation, and account for the observed density of lapetus provided that this 
satellite formed in situ (Mosqueira and Estrada, 2005). Although we focus on lapetus, our 
model applies to the origin of the regular satellites in general^. The model we present here 
fits with the gaseous SEMM satellite formation model (Mosqueira and Estrada 2003a,b; 
hereafter MEa,b; see also Estrada et al. 2009 and Mosqueira et al. 2009). This framework 
is consistent with the core accretion planet formation model in which Jupiter and Saturn 
form ~ 3 - 5 Myr after CAIs (Hubickyj et al. 2005; Dodson- Robinson et al. 2008). Note 
that in alternative gas-poor or gas-starved satellite formation models planetesimals do not 
traverse enough subnebula gas to ablate significantly. Thus, mass delivery by ablation is not 
available to these models whether they consider planetesimals as a source of solids (Estrada 
and Mosqueira 2006) or not (Canup and Ward 2002). 

For rock/water-ice fractionation to take place, and our lapetus' model to obtain, the 
first generation of planetesimals in the Jupiter-Saturn region must incorporate significant 

■^For instance, the puzzling observation that primordial ^^Ar is enhanced in the atmosphere of Jupiter 
(Atreya et al. 1999), but only present in trace amounts in the atmosphere of Titan, as reported by the 
GCMS instrument aboard Cassini (Niemann ct al. 2005), may fit within a framework of planetesimal 
ablation followed by subnebula re-condensation of some volatiles but not others, but we leave this subject 
for further work. 
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quantities of ^^Al, which depends on their time of formation. The closest analog to the 
planetesimals forming in the Jupiter-Saturn region is the asteroid belt. The ages of asteroids 
indicate that they formed with a spread of 1 — 3 Myr after Ca-Al-rich inclusions (CAIs) 
(Scott and Krot 2005). The asteroid belt may have been cleared largely by scattering due 
to planetary embryos and a migrating Jupiter (Bottkc ct al. 2005). Yet, the fossil evidence 
provided by asteroids may not be directly linked to the first generation of planetesimals. 
(It is worth noting that the parent bodies of iron meteorites may have formed 1 — 2 Myr 
earher than those of the ordinary chondrites; Baker et al. 2005.) Instead, main belt asteroids 
arc the product of a few My of accretion before most of the mass in the belt was depleted 
(Weidcnschilling 2009). Indeed, as runaway growth proceeds to completion most small bodies 
would be ground down by embryos. The degree of coUisional evolution in asteroidal objects 
depends on the size of the planetesimals, which may be used to provide a constraint of ~ 
10 — 100 km on primordial planetesimal sizes (WeidenschiUing 2009). Thus, neither asteroids 
nor meteorites provide a direct sample of this population. Nevertheless, iron meteorites may 
have come from ~ 20 — 200 km differentiated parent bodies that underwent catastrophic 
collisions (Chabot and Haack 2006). 

The timescale for the formation of planetesimals depends on the nebular environment 
early on. Assuming efficient sticking, the timescale for the formation of planetesimals is 
very short indeed r ~ 2000 (a/lAU)^/^ years, where a is the radial distance from the Sun 
(WeidenschiUing 2000). CoUisional growth appears to be efficient for small particle sizes and 
impact speeds (see review by Dominik et al. 2007); yet, growth past the decoupling size ~ 1 
m in the presence of turbulence is questionable (Mosqueira 2004; Cuzzi and Weideschilling 
2006). The problem is especially severe in light of the short clearing times of ~ 1 m size 
objects. It can be argued that in the presence of turbulence planetesimals can bypass this 
size range entirely (Johansen et al. 2007; Cuzzi et al. 2008), however there is no agreement 
regarding the specific mechanism to accomplish this leap - the simulations of Johansen et al. 
(2007) assume a particle size that is already partially decoupled from the gas, whereas the 
simulations of Cuzzi et al. (2008) do not yet show that a large planetesimal > 10 km can 
actually form (even sporadically) from chondrule-like ~ 1 mm sized particles. While a recent 
study concludes that asteroids are born big (Morbidelli et al. 2009), given uncertainties in the 
strength of planetesimals and the role of coUisional debris, the asteroid size distribution, and 
the timing of the excitation caused by the migration of Jupiter, planetesimal sizes 10 — 100 
km appear consistent with the data. In this regard, the possibility should be explored that 
smaller planetesimals are disrupted by planetary embryos prior to the final stages of asteroid 
accretion. 

The alternative that quiescent regions can help to form planetesimals has been considered 
in a number of publications (Krctkc and Lin 2007; Brauer et al. 2008; Lyra et al. 2009). 
The stumbling block here may be an embarrassment of riches, i.e., in a laminar disk it is 
too easy to make planetesimals. Consequently, dust might coagulate into larger objects too 
fast to explain the observations of protoplanetary disks indicating the presence of small < 3 
//m grains (e.g., Dullemond and Dominik 2005; Cuzzi ct al. 2008; Birnstiel et al. 2009), 
the evidence for delayed (~ 1 Myr) grain growth based on observations of SEDs of young 
disks (Beckwith et al. 2000; Currie et al. 2009), or the lack of pervasive differentiation of 
main belt asteroids (Cuzzi et al. 2008). However, dust may be replenished by weak inflow 
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onto the protoplanetary disk (Dominik and DuUemond 2008), by viscous mixing from out- 
side the "dead-zone" (Birnstiel et aL 2008), or by coUisional grinding; in addition, main 
belt asteroids do not constitute a fossil record of the first generation of planetesimals. The 
degree of activity in the "dead-zone" remains to be clarified (Turner et al. 2007; Bai and 
Goodman 2009), but there is general agreement in the literature that disk quiescence facili- 
tates planetesimal formation (e.g., Youdin and Shu 2002; Cuzzi and Weidenschilling 2006). 
Since the protoplanetary disk may remain turbulent and variable during its initial stages (the 
"FU-Ori" epoch) lasting a few 10^ years (Bell et al. 2000), we take < 10^ years to be the 
planetesimal formation timescale, implying that these planetary building blocks incorporate 
significant quantities of short-lived radioactive elements. 

In Sec. 2 we discuss the framework and model parameters for the computation of the mass 
delivery to the circumplanetary disk. In Sec. 3 we describe the compositional constraints that 
our model seeks to address, and the physical processes used to that end. We briefly discuss 
the thermal evolution of planetesimals in the Jupiter-Saturn region due to the presence of 
SLRIs in Sec. 4. We close with our conclusions in Sec. 5. 

2 Model parameters 

Following giant planet formation planetesimals in its feeding zone undergo collisional grind- 
ing. In the Jupiter- Saturn region the collisional timescale for kilometer-sized objects is 
similar to the ejection timescale ~ 10^ yrs, so that a fraction of the mass of solids will be 
fragmented into objects smaller than 1 km (Stern and Weissman 2001; Charnoz and Mor- 
bideUi 2003). The collisional cascade facilitates planetesimal dehvery to the circumplanetary 
disk because smaller planetesimals are easier to capture. An estimate of the size that will 
be fragmented may be obtained by equating the ejection time Tgjerf ~ 0.1Q~^[Mq/Mp)'^, 
where fl is the angular velocity, Mp is the planet's mass and Mq is the Sun's mass, to the 
collision time ~ pgr/flag, where Ps is the density of the solids, r is the particle size and 
as is the surface density of the solids. This yields particles as large as 10 km at Jupiter (and 
larger for Saturn) for surface densities a few times the minimum mass (Hayashi 1981). 

Switching to the subnebula, in the SEMM model (MEa,b) one expects a factor of ~ 10 
enhancement in solids over cosmic mixtures, resulting in a gas surface density of ~ 10^ g cm~^ 
(see Sec. 2.2), which is consistent with and quantitatively constrained by the gap-opening 
condition for Ganymede and Titan in a inviscid disk with aspect ratio ~ 0.1 (Rafikov 2002; 
MEb), and the Type I migration of full-sized satellites in such a disk (MEb; Bate et al. 
2003). If we use a gas surface density of 10^ — 10^ g cm~^ for the Jovian and Saturnian 
subnebulae, then a planetesimal of density ""^ 1 g cm~^ will encounter a gas column equal 
to its mass if its radius is in the range 0.1 — 1 km, resulting in ablation and delivery to the 
circumplanetary disk. 

Our aim here is to provide a link via ablation of planetesimal fragments in the Jupiter- 
Saturn dynamically active region to the circumplanetary disks of these two planets, seeking 
to account for the mass and angular momentum budgets and compositions of the regular 
satellites. Tying solar nebula planetesimals to subnebula satellitesimals involves three dif- 
ferent aspects: first, characterizing the properties of the first generation of planetesimals in 
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terms of sizes and degree of heterogeneity; second, quantifying the coUisional evolution of 
the planetesimal swarm following giant planet formation; and third, delivering planetesimal 
fragments to the circumplanetary disks of Jupiter and Saturn. 

2.1 Planetesimal fragmentation and the initial size distribution 

The strength of planetesimals is an unknown quantity. A strain-rate scaling law by Housen et 
al. (1991) suggests that kilometer-sized objects are the weakest. Davis and Farinella (1997) 
find a constant S — 3 x 10^ erg cm~^, where S is the impact strength (i.e., the energy per 
unit volume for shattering fifty percent by mass of the parent body) for crushed icy bodies, 
in agreement with Ryan et al. (1999). Benz and Asphaug (1999) use a 3D Hydrocode to 
simulate 3 km/s impacts and find that 100 m objects may be weakest, with kilometer sized 
planetesimals strengthened by gravity. How much planetesimal mass fragments following 
the coUisional cascade and is delivered depends on the planetesimal strength (Charnoz and 
MorbideUi 2003). 

The initial size distribution is also unknown. If run-away growth takes place the resulting 
size distribution may be steep. Wcthcrill and Stewart (1993) find a bimodal size distribution 
with a fragmentation tail of objects less than 1 km with power law exponent q ~ 3.5 and size 
distribution for objects larger than that with ~ 5.5, where the differential size distribution 
{dN/dr oc r~^) approaches a value of g ~ 3.5 in the case of a size independent fragmentation 
model (e.g., Dohnanyi 1969). In his coagulation simulations Weidenschilling (1997) finds 
that most of the mass resides in objects ~ 10 km in a timescale of 10^ — 10^ years in the 
outer nebula. Most studies of planet formation start with most of the mass of solids in 
objects of size ~ 10 km, as suggested by numerical simulations (e.g., Kenyon and Luu 1999). 

Unless otherwise stated, we consider values oiq — 3.5. We set the lower size cut-off at 1 m. 
The reason is that smaller objects are protected from further coUisional grinding by gas drag, 
so that below this size collisions arc once again accretive as long as the nebula turbulence is 
weak. We vary the upper size cut-off between 10 — 100 km in lieu of more detailed coUisional 
fragmentation simulations, which we leave for further work. For simplicity, we also assume a 
fully fractionated population of icy and rocky (which includes iron) planetesimal fragments, 
each of which has the same particle size spectrum (see Sec. 4) . 

2.2 Disk properties 

Gas flowing into its Hill sphere forms a disk around the protoplanet. The timescale for 
Jupiter and Saturn to clear a gap in between ~ 10"^ years sets the time for the end of gas 
accretion. The timescale for envelope collapse is given by the Kelvin- Helmholtz time of 
10^ — 10^ years (Hubickyj et al. 2005), following which a circumplanetary disk forms. Prior 
to gap-opening the giant planet accretes gas with semi- major axis originating from ~ Rh of 
its location, where = a(Mp/3MQY^^ is the planet's Hill radius, centrifugal balance yields 
a characteristic disk size of ~ Rh/50 (Stevenson et al. 1986; MEa). For Jupiter and Saturn 
these radii are located close to the positions of Ganymede and Titan. After gap-opening 
accretion continues through the planetary Lagrange points, and the estimated characteristic 
disk size is larger ~ Rh/^ (MEa; Estrada et al. 2009; Ayliffe and Bate 2009), which we take 
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to be the radial extent of the circumplanetary disk. This radial size corresponds to ~ 150 
Rj, and ~ 200 R5, where Rj and Rs are the radii of Jupiter and Saturn, respectively. 

MEa divide the circumplanetary disk into inner and outer regions. For Jupiter, we 
compute the solids-enhanced minimum mass (SEMM) gas densities in the inner and outer 
disks based on the solid mass required to form lo, Europa (both re-constituted for missing 
volatiles) and Ganymede in the inner disk, and Callisto in the outer disk. Inside of the 
centrifugal radius the surface gas density is S ~ 10^ — 10^ g cm~^, which corresponds to 
pressures ~ 0.1 bar. Outside the centrifugal radius Rc ~ 30 Rj, the gas surface densities 
are in the range 10^ — 10^ g cm~^. We apply the same procedure to the circumplanetary 
disk of Saturn by employing the masses of Titan and lapetus to set the inner and outer disk 
masses, respectively. The transition region has a width > 2Hc, where He is the subnebula 
scale-height at the centrifugal radius. This choice ensures that the gradient in gas density is 
not so steep as to lead to a Rayleigh- Taylor instabihty (e.g., Lin and Papaloizou 1993). 

The subnebula gas surface density profile is taken to be 

{'^a(Ra/R), R < Ra] 

E„(i?„/it:) - Ra<R<R,; (1) 

e 2ri +1 
^b{Rb/ R^: R > Rb- 

Once a disk forms, the cooling timescale depends on the opacity. As particles grow, the 
opacity of the nebula decreases, allowing the gas to cool and accretion of ice-rich satellites to 
take place. We connect the disk's temperature profile to the planetary luminosity at the tail 
end of giant planet formation (Hubickyj et al. 2005). We use a heuristic temperature profile 
of the form T = 3600 Rj/i? for Jupiter (e.g., Lunine and Stevenson 1982) and T = 2000 Rs/-R 
for Saturn. The outer disks of Jupiter and Saturn have roughly constant temperatures in 
the range of 70 — 130 K for Jupiter and 40 — 90 K for Saturn, depending on solar nebula 
parameters. 

3 Regular Satellite Mass and Composition Constraints 

The compositional gradient of the Galilean satellites may provide a probe to the environ- 
ment in which they formed (Estrada et al. 2009 and references therein). Here we focus on 
the large, outer regular satellites of each satellite system: Ganymede and Callisto in the 
case of Jupiter; and Titan and lapetus for Saturn. Planetesimal break-up in tandem with 
delivery via ablation of planetesimal fragments crossing the subdisk provides a framework 
for understanding the mass budget and compositions of regular satellites compared to that 
of solar nebula planetesimals. 

3.1 Thermal Ablation of Disk Crossers 

Thermal ablation occurs because friction of the fast moving heliocentric interloper as it 
crosses the circumplanetary gas disk heats up the body. At low gas densities impactors lose 
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mass and energy through ablation. For large kilometer-sized bolides mechanical destruction 
rather than thermal ablation may dictate the fate of the object (e.g., Zahnle and Mac 
Low 1994). Ablation of planetesimal fragments < 1 km may deliver the bulk of the solids 
needed to form the satellites. For objects in that size range, the heat transfer coefficient 
Ch ~ 0.1 may be reasonably obtained from observations of terrestrial meteorites (Bronshten 
1983). The rate at which energy is transferred to the planetesimal per unit area is given by 
E ~ Q.hCHpv^-i where p is the gas density, v is the speed of the planetesimal through the gas, 
and some planetesimal flattening (which increases its cross section) takes place (e.g., Zahnle 
1992; Chyba et al. 1993). Ablation thus results in delivery of material to the circumplanetary 
disk. Gas drag can also result in capture of material. 

Ignoring conduction into the interior and ablation, we can obtain the surface temperature 
Tsurf by balancing this heating and radiative cooling easBTg^j.f, 

Tsurf-^Ji^pv' + TA''\ (2) 



V 8(755 ' 

where e = 0.5 is the emissivity, Tq is the subnebula background temperature, asB is the 
Stefan-Boltzmann constant, v is the speed of the planetesimal through the gas disk, and 
/ is a parameter that measures the degree of flattening. It is important to keep in mind 
that this quantity reflects the energy regime not the actual surface temperature unless it is 
low enough that ablation can be ignored. Using parameters appropriate at the location of 
lapetus ~ 5 X 10^ cm (a cold population at infinity), / = 4 and p = T,/2H ~ 7 x 10"-^° 
g cm~^, where E is the gas surface density and H — c/Q is the scale-height, c the sound 
speed and Q = a/ GMp / is the angular velocity in the circumplanetary disk, we obtain a 
surface temperature of Tgurf ~ 900° K, which is sufficient to melt and vaporize icy objects 
but not rocky objects. A value of / = 1, which is suitable to the outer disk due to the low 
ram pressures there, yields a temperature of ~ 600° K. In contrast, at Titan (^; ~ 8 x 10^ 
cm p ~ 3 X 10"''' g cm~^) and Callisto (■?; ~ 1.1 x 10^ cm p ~ 2 x 10"''' g cm~^) the 
surface temperature can exceed 4000° K (for / = 4), which is enough to melt and vaporize 
rocky objects as well. Hence, the temperature of objects crossing the disk places lapetus in 
a separate regime from other regular satellites, i.e., at its location in the disk rocky objects 
do not reach temperatures high enough to ablate. 

We can estimate the change in radius of an icy planetesimal due to either melting or 
vaporization as it crosses the gas disk in a time At ~ 2H/v (ignoring gas drag for now) 
(Podolak et al. 1988) 



\dt / yap,melt ^Pv^Oc-m 



where pp is the planetesimal density, T^. = 648° K is the critical temperature for vaporization, 
and = 273° K is the melting temperature, i^om = 6.0 x 10^ erg g~^ is the energy required 
to raise the temperature of 1 g of ice from the background temperature of the nebula ~ 100° 
K (where it originated) to plus the latent heat of melting Hf, and Eqc — 2.8 x 10^° erg 
g~^ is the energy needed to raise the temperature of 1 g of icy material for 100° K to Tc (so 
that vaporization becomes energy limited) plus the latent heat of vaporization Hy. We find 
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that the amount of ice melted or vaporized at lapetus {At ~ 3 x 10^ s) is in the order of 
meters (for / = 1), whereas at Titan and CaUisto {At ~ 4 x 10^ s) it is about a kilometer 
(for / = 4). It should be noted that at lapetus (E ~ 100 g cm~^) meter-sized objects (with 
density ~ 1 g cm~^) traverse a column of gas equal to their mass, whereas at Titan and 
Callisto (S ~ 10'^ g cm^^) the same is true for kilometer-sized objects. Hence, ablation 
of planetesimal fragments may provide the solids needed for the formation of the regular 
satellites provided that a significant fraction of the planetesimal mass resides in fragments 
with sizes in the kilometer to meter size range. 

One can ask whether a particle will ablate versus breakup when crossing the disk. At 
Titan, the dynamic pressure is Pdyn = ~ 10^ dynes cm~^. At lapetus, we find 

Pdyn ~ 100 dynes cm^^. Therefore, if a meter-sized or larger particle with strength S ~ 10^ 
dynes cm~^ crosses the disk icy material is left behind by ablation (at Titan and lapetus) 
or by gas drag capture (at Titan), but probably not breakup (at least at lapetus). It is also 
important to consider the sensitivity of our results to the flattening parameter /. In this 
regard, it should be noted that the overall dependence of ablation on / is not too strong, 
and that at the location of lapetus the ram pressure is low enough that / = 1 is a good 
approximation. 

To find out the mass budget and composition of the subnebula resulting from the ablation 
of an unmixed population of rocky and icy disk crossers, we have to tackle the thorny issue of 
whether the rate of ablation in Eq. 3 is dominated by melting or vaporization. It is probably 
fair to say that most models of the ablation process focus on vaporization (e.g., Moses 1992). 
However, at lower temperatures the rate of ablation may be affected by the viscosity of the 
melt and other factors (see Podolak et al. 1988 for a discussion). Yet, a detailed model 
incorporating the full complexity of the problem for a broad range of parameter space is 
premature in that it can obscure our conclusions. In fact, as we show below, our main results 
are robust and not dependent on the specifics of the ablation mechanism. For this reason, 
we leave development of a more detailed ablation model for later work, and we employ 
the classic approach of Bronshtcn with minor modifications. We calculate the amount of 
material ablated by solving coupled equations of motion and ablation (Bronshten 1983) for 
a planetesimal of mass m which is assumed for simplicity to cross the circumplanetary disk 
perpendicular to the disk plane: 

= -lcDA{m)p{X, Zy - mn''{X)z; (4) 
at 2 



dm 
~dt 



C 



H 



2E, 



Oc.m 



-A{m) 



'p{X,Zy-^-^{eTl^-T^) 



C 



H 



for Tg^j-f > Tc^jji, 



^ _ n f T < T 

'J lOr J- surf — -'- c,m-i 



(5) 



where T^^r/ is obtained from Eq. 2 with / = 1, C^j = 0.44 is the drag coefficient, A is the 
planetesimal cross-section, X — R/Rp is the distance from the planet in units of planetary 
radii Rp, Z — z/V^H, p{X, Z) — po{X)e~^ is the gas density of the subnebula. Here we 
set the flattening factor f — 1 because we are particularly interested in the outer disk, where 
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the ram pressure is low. Note that for Tsurf < T'c,m the ablation equation goes to zero even 
though the planetesimal can still be captured by gas drag. That is, regions of the disk such 
that the disk crosser never reaches the ablation temperature Tc^m may only be delivered to 
the circumplanetary disk by gas drag capture. Fits to observational data of meteoroids in 
the atmosphere yield values of Ch — 0.1 for low values of the effective heat of ablation ~ 10^° 
erg g~^ (Svetsov et al. 1995), which leads us to use £^om = 3-0 x 10^° erg g~^ and an ablation 
temperature of T^^ — 1800° K for rock. For ice we consider both melting (-Eom — 6.0 x 10^ 
erg g-i; = 273° K) and vaporization {Eqc = 2.8 x 10^° erg g-^; = 648° K) while 
keeping Ch — 0.1 constant. Wc stress that these parameter choices are quite conservative 
and intended to test the sensitivity of our results, as significant water-ice sublimation will 
take place for Tsurf < Tc,m- Thus, a more realistic treatment is likely to deliver significantly 
more water to the outer disk. 

We solve Equations 4 and 5 using a fourth order Runge-Kutta scheme for coupled equa- 
tions. In Fig. 1 we plot final velocity relative to initial v final /vinit-i final planetesimal size 
1^ final, distance travelled through the disk z/H, and delivered mass fraction due to ablation 
and gas drag capture as a function of the distance from Jupiter, X, for different particle sizes. 
Here the gas surface density corresponds to that of a SEMM model = 3x 10^/Xg cm~^ 

(without a transition). These results use the more conservative ice vaporization values for 
the ablation temperature Tc and heat of ablation Eqc. If the planetesimal does not com- 
pletely ablate, it may successfully pass through the disk. However, planctcsimals may fail to 
traverse the disk because they are ablated before they do, or because of gas drag capture. If 
particles are ablated down to a size of 1 cm, they are taken to have totally ablated. Close 
to the planet even ~ 1 km planetesimals (bold solid lines) are completely ablated before 
reaching a final velocity v final — (top left panel). Farther out in the disk, a leftover plan- 
etesimal of size r final (top right panel) emerges from the disk after crossing a distance of 
2H (bottom left panel). As expected, the mass delivered to the disk by either capture or 
ablation decreases with distance from the planet (bottom right panel). 

In the top left curve of Fig. 1 we also plot the escape velocity Vesc = \/ 2GMp / (RpX) 
as function of distance (light solid line) for reference. It is important to point out that here 
(and in the rest of this paper) our capture condition is defined locally. That is, we only 
consider that the particle has been locally captured if v final = 0. A planetesimal can emerge 
from the disk with a speed too slow to escape the gravitational potential well of the planet, 
so that it will be delivered either to the circumplanetary disk or the planet. For instance, 
10 m particles (short dashed lines) have final speeds that are less than the escape speed out 
to the edge of the disk. Therefore, these particles will be delivered to the circumplanetary 
disk, but not always locally. Here we require v final — for capture since in this case the 
planetesimal mass is delivered at the location of passage. This is conservative in the sense 
that it underestimates the mass delivered to the circumplanetary disk. But it is likely to 
provide a better estimate of the ice/rock ratio as a function of position because particles with 
significant non-zero exit velocities can potentially be delivered anywhere inside the crossing 
radius. We leave a more realistic treatment of multiple passes through the disk for later 
work. 

In Fig. 2 we do the corresponding calculations for the case of rock. Note that in the case of 
rock we employ the rock melting values for the ablation temperature Tm and heat of ablation 
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Figure 1: Plot final velocity relative to initial v final /vinit, final planetesimal size r final, distance 
travelled through the disk z/H, and delivered mass fraction due to ablation and gas drag capture 
as a function of the distance from Jupiter, X, for several initial planetesimal sizes (indicated by 
the solid and dashed curves). The gas surface density is taken to be that of the SEMM model, 
= 3 X 10^ /X g cm^^ (without a transition). We take the initial velocity as the interloper 
crosses the disk to be vo ~ GO/s/X km s"^. The planet's escape velocity as a function of distance 
Vesc is plotted in the top left panel. 

Eonj,. For rock the mass delivered for a given particle size (bottom right panel) drops faster 
than for ice despite our conservative choice of parameters for ice ablation. Consequently, 
the composition of disk is expected to be enriched in ice. Rock mass delivery in the outer 
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Figure 2: Plot final velocity relative to initial Vfinai/vinit, final planetesimal size r final, distance 
travelled through the disk z/H, and delivered mass fraction due to ablation and gas drag capture 
as a function of the distance from Jupiter, X, for several initial planetesimal sizes (indicated by 
the solid and dashed curves). The gas surface density is taken to be that of the SEMM model, 
= 3 X IQ^/X g cm~^ (without a transition). The planet's escape velocity as a function of 
distance Vgsc is plotted in the top left panel. 

disk can be dominated by capture. This can be seen by the behavior of 1 m particles (long 

dashed lines): gas drag slows down these particles sufficiently fast that almost no ablation 
takes place beyond X = 60 so that r final = 1 m (top right panel). Yet gas drag is effective 
out to X = 70 despite the low surface density of the outer disk so that v final = (top 
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Figure 3: Number of particles with at least one crossing inside a given distance from Jupiter. 
Simulations involved 5000 particles with total mass of 10 Mq. We estimate that roughly ~ 7 
pass within 0.2 Rh- 

left panel) and the full particle mass is delivered by gas drag capture out to this location 
(bottom right panel). Outside this location essentially no rock is delivered. Because ice is 
still ablated and captured outside this location (Fig. 1), one can expect the cold outer disk 
to be water-ice enriched following re-condensation. 

We follow the dynamical evolution of a 10 mass disk (where is an Earth mass) 
of 5000 particles spread between 4 and 8 AU starting on low eccentricity and inchnation 
orbits. Jupiter is placed at 5 AU and Saturn is at 7.5 AU (packed configuration), both 
with low values of eccentricity and inclination. In Figures 3 and 4 wc plot the number of 
particles that crossed at least once within a given distance of Jupiter and Saturn in 10^ 
years. Not surprisingly, a large fraction of the total number of particles in the simulation 
(5000) cross at least once within 0.1 of the Hill sphere of Jupiter or Saturn in the first 
~ 10^ yrs. This is because particles are scattered during close encounters with the giant 
planets. Several caveats apply to these results. For one, collisional fragments, which are 
easier to deliver to the circumplanetary disk, may collide frequently, which would change 
their dynamics. This process is not modeled in the results discussed above, and is left for 
future work. But we stress that collisions cannot prevent collisional debris from crossing 
the circumplanetary disk of either giant planet - particles of any size at the gap edge and 
beyond are able to penetrate the Hill sphere of the giant planet during close encounters. In 
some regimes of interest, such as cases in which nebula gas drag and collisional damping 
is included and metcr-sizcd particles play an important role, one needs to treat the full 
fragmentation-coagulation problem. However, it is presently unclear how much mass the tail 
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Figure 4: Number of particles with at least one crossing inside a given distance from Saturn. 
Simulations involved 5000 particles with total mass of 10 M^. We estimate that roughly ~ 5 
pass within 0.2 Rh- 



end of the mass-spectrum contains. 

We calculate the total captured mass of the ice and rock components Mj{X) dehvered 
within a radius X — R/Rp to the circumplanetary disk in time r using (e.g., Estrada and 
Mosqueira 2006) 



(X) = 2t:tRI Fj {X')X' dX'- 
Fj{X) = voil + ie/X) / Cj<^j{X,m)m^-Pdm 

«/ TTlrnAn. . n 



(6) 



where Fj is mass inflow deposition rate per unit area of component j, Vq is the speed at 
infinity, p — (1/3) (g' + 2), rrimin,] ^-nd rrimax^j ^-re the lower and upper bounds of the mass 
distribution of component j, = GMp/Rpv"^ is the Safronov parameter, / ~ 1 — 2 is a 
constant, and $j(X, m) is the fraction of a planetesimal with mass m that is delivered to 
the disk at X, as obtained from Eqs. 4 and 5. The Cj are normalization coefficients given by 



{2-p)Mpi 
27rRpVoT 
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for ice; 
for rock ' 
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where Mpi is the mass of planetesimals passing through the disk, Xr is the mass fraction of 
rock, and for our disk size we take — 0.2Rh/Rp. The surface density of components j, 
Ej, sohds-to-gas-ratio, and the ice mass fraction as a function of radius are obtained from 
the Fj{X). 

Pending more detailed scattering/fragmentation simulations, we assume a power- law 
mass distribution with slope q = 3.5 and an upper and lower particle size cut-offs of 
rmax = {'^rrimax/^T^PpY^^ = 200 km and Tmin = {'^rrirain/^T^PpY^^ = 0.001 km for both rock 
and ice. We allow for a differentiated population of icy and rocky fragments with a 30/70 
ratio by mass. Based on our A^-body simulations, we estimate that about Mpi ~ 7 cross 
at least once within 0.2 Kh- It should be noted that particles ~ 1 km that cross the disk 
at > 0.1 Kh are not in general ablated so they may cross again within 0.2 Kh before being 
ejected. We take the Safronov parameter B = 25 representative of a fairly cold population. 
The factor I accounts for the different geometry of a flux of impactors through the subnebula 
rather than the planet itself (MorfiU et al. 1983; Cuzzi and Estrada 1998; Charnoz et al. 
2009). The value for / should be obtained from numerical simulations. For now wc use / = 2, 
but the results are not sensitive to this parameter. The background nebula temperature is 
set to 130 K. For our subnebula parameters, we choose = 2 x 10*^ g cm~^, = 5 x 10^ 
g cm~ ^, Rg = 1 5 Rj, Rb = 25 Rj, and n — 1. The velocity of the planetesimal is given by 

In Fig. 5 we show the mass delivery by ablation and gas drag capture of planetesimal 
fragments crossing the SEMM circumplanetary disk. The top left panel indicates the size 
of the largest particle that will be delivered to the disk in full by a combination of ablation 
and capture of the leftover object (dotted curves). It can be seen that a higher proportion 
of mass is delivered by capture for rock than for ice (especially for the ice melt curve). This 
is because less ablation makes the rock particle cross section and gas drag force larger for 
a given size. The top right panel corresponds to the solids/gas fraction that results from a 
cold population of planetesimal fragments as a function of position in the disk. The surface 
gas density is also indicated. Both the dotted and the dashed curves exhibit a bump during 
the inner to outer disk transition. This jump in the concentration of solids occurs because in 
the transition region the surface gas density of the disk decreases faster than the decrease in 
the mass delivery. Consequently, the concentration of solids increases. The subsequent outer 
disk rise in the dotted curve and bump in the dashed curve are related to the l/R"^ surface 
density profile in the outer disk (compared to a 1/i? profile in the inner disk) together with the 
transition to an isothermal temperature profile in the outer disk (compared to a 1 /i? profile 
in the inner disk) . The dashed curve drops towards the outer edge because the temperature 
of the planetesimal becomes too low for significant vaporization either of rock or ice (whereas 
the temperature is still high enough to melt water-ice). The bump in the concentration of 
solids at the transition between the inner and outer disk is tantalizing in that it can help to 
accrete satellites close to the present location of Ganymede (and Titan for Saturn's disk). 
But it is premature to say that this feature is significant. The bottom left panel indicates 
the resulting ice/rock fraction of the disk. The bottom right panel indicates the total mass 
delivered within a radial distance R. The region between the bars contain enough mass to 
make each of the indicated satellites (reconstituted in the case of lo and Europa). The mass 
delivered was obtained by using the ice melt results for Tsurf < Tc and the vaporization 
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results for larger Tgurf- These results indicate that enough mass can be delivered to account 
for the reconstituted Gahlean satellites. But we stress that the mass delivered to the disk 
is dependent on the parameters chosen and the capture condition v final — 0. On the other 
hand, the capture condition is a conservative choice in that significantly more mass will be 
delivered to the circumplanetary disk than we assume here. We reserve a parameter study 
for later work. 

The two ice curves shown correspond to ablation by melting and vaporization in Eq. 5. 
As one would expect, more ice is delivered when ablation turns on at the lower tempera- 
ture of melting versus vaporization Tg. The good agreement with the mass budget and 
compositions of Ganymede and Callisto would appear to argue in favor of the vaporization 
curves. However, this conclusion is premature for several reasons. First, these calculations 
make the assumption that the population of disk crossers is completely fractionated, which 
is unrealistic even if primordial planetesimals were fully differentiated (see Sec. 4), as some 
mixing would be expected during coUisional disruption. A more realistic treatment would 
lower the ice/rock ratio of the entire disk. Two, vaporization of water-ice actually starts 
below the melting temperature Tm = 273° K (whereas for rock temperatures 200° K above 
melting Tm = 1800° K are needed before significant vaporization takes place; Podolak et al. 
1988). Third, the total mass dehvered to the circumplanetary disk is dependent on parame- 
ter choices. However, the result that the outer disk is icier than the inner disk is independent 
of the specifics of the ablation model. In fact, given a fully fractionated population crossers, 
the outer disk receives material mostly by ablation of ice, and essentially no rock either by 
ablation or capture. Also, these results lead one to expect that the circumplanetary disk is 
enriched in ice over solar mixtures for a broad range of parameter choices. 

Similarly for Saturn, we assume a power-law mass distribution with slope q — 3.5 and 
allow for a differentiated population of icy and rocky fragments with a 30/70 ratio by mass. 
Based on our simulations, we estimate that about Mpi ~ 5 cross at least once within 
0.2 Kh of Saturn. We let = 100 and an upper and lower particle size-cutoffs of Vmax = 50 
km and Tmin = 0.001 km for both rock and ice based on the larger fraction of rubble and 
lower velocity dispersion expected for Saturn. The background nebula temperature is set 
to 90 K. For our subnebula parameters, we chose = 2 x 10^ g cm~^, E^, = 2 x 10^ g 
cm~^, Ra = 20 R5, Rh = 58 R5, and n = 3. As we can see from Fig. 6, rocky fragments 
in the 0.001 — 1 km are not either captured or ablated, whereas icy fragments are delivered 
mostly by ablation. Thus, delivery of solids via ablation of a non-homogeneous population 
of planetesimal fragments leads to an enrichment of water-ice in Saturn's outer disk. 

In summary, these results show that it is possible to dehver enough mass to account for the 
regular satellites of Jupiter and Saturn. Furthermore, ablation and the capture cross-section 
of meter-sized objects can lead to ice/rock fractionation and account for the composition 
of lapetus, as well as those of Ganymede, Callisto and Titan. This is the only explanation 
currently available for lapetus' icy composition, and may also indicate that regular satellites 
overall may be depleted in rock with respect to solar composition mixtures, which might 
explain why their densities are lower than that of similar size Kuiper belt objects. We stress 
that while the amount of mass delivered to the circumplanetary disks of Jupiter and Saturn 
is dependent on parameter choices, such as the upper size cut-off and the slope of the mass 
spectrum, the velocity dispersion at infinity and the capture condition, we found that the 
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Figure 5: Left upper panel: The planetesimal size that gets ablated (solid and dashed curves) and 
captured (dotted curves) as function of radial distance. The dotted horizontal line corresponds to 
the lower size cut-off. Right upper panel: The solids to gas ratio for the two ice models (dotted 
and dashed curves) as a function of radial distance from Jupiter. The bold solid curve indicates 
the gas surface density of the disk in g cm^^ (right side). Left lower panel: The ice to rock ratio 
as a function of radial distance. The solid dots correspond to Ganymede (inner disk) and Callisto 
(outer disk). The ice mass ratio of solar nebula planetesimal fragments crossing the disk is taken 
to be 0.3. Lower right panel: Plot of the mass delivered as a function of distance from Jupiter. 
The two models correspond to melting and vaporizing ice. 
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Figure 6: Left upper panel: The planetesimal size that gets ablated (solid and dashed curves) and 
captured (dotted curves) as function of radial distance. The dotted horizontal line corresponds to 
the lower size cut-off. Right upper panel: The solids to gas ratio for the two ice models (dotted 
and dashed curves) as a function of radial distance from Saturn. The bold solid curve indicates the 
gas surface density of the disk in g cm~^ (right side). Left lower panel: The ice to rock ratio as 
a function of radial distance. The solid dots correspond to Titan (inner disk) and lapetus (outer 
disk). The ice mass ratio of solar nebula planetesimal fragments crossing the disk is taken to be 
0.3. Lower right panel: Plot of the mass delivered as a function of distance from Saturn. The two 
models correspond to melting and vaporizing ice. 
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efficiency of delivery is high enough that the mass budget of the satelhtes can easily be 
accommodated. We find that the water-enrichment of the outer disks of Jupiter and Saturn 
is a robust feature given a differentiated population of icy and rocky interloper fragments. 

While full ice/rock fractionation is unlikely, partial separation can account for lapetus' icy 
composition. Such a population results from the fragmentation of a partially differentiated, 
^^Al-heated ffist generation of ~ 10 km planetesimals (Mosqueira et al. 2009). Indeed, it is 
likely that most of the mass in the ffist generation of planetesimals resided in objects 10 — 100 
km in the first 10^ — 10^ years, so that these objects incorporated significant amounts of ^^Al, 
leading to partial differentiation, as we show below. 



Table 1: Physical Parameters 



Parameter 


Ice 


Rock 


fx (g mole~^) 


18 


50 


Ps (g cm"^) 


0.92 


3.4 


Eoc (erg g"^) 


2.8 X 10^° 


8.08 X 10^° 


Eom (erg g-^) 


5.94 X lO'^ 


3.04 X 10^° 


Hf (erg g-^) 


3.33 X 10*^ 


1.43 X 10^° 


Tm (°K) 


273 


1800 


Tc ("K) 


648 


4000 



4 Differentiation of primordial planetesimals 

The thermal evolution of comets (Prialnik et al. 2008) and asteroids (Ghosh et al. 2006) 
powered by short-lived radionuclides has been treated in an extensive number of publications, 
and it is not possible for us to do it justice in this brief section. Our aim here is to get a 
handle on the degree of differentiation of ~ 10 km primordial planetesimals forming in the 
Jupiter-Saturn region. This population of objects would have formed in a faster timescale 
and in a hotter region of the nebula than their Kuiper belt counterparts, yet still outside 
the snowhne. With the possible exception of ~ 100 km Phoebe, which may itself be at least 
partially differentiated (Johnson et al. 2009), and whose composition appears to match that 
of other outer solar nebula objects (Johnson and Lunine 2005), there are no known examples 
of this population. 

Our simplified treatment follows along the same lines as that of other works only applied 
to our region of interest. We leave a more accurate study for later work (in collaboration 
with M. Podolak). We solve the energy balance equation for a body of size r up to the 
melting temperature in a time tm (e.g., Jewitt et al. 2007) 

.rr^^dT 3k(T)T 
c,{T) ^ = E ^^^o,jHj - (8) 

where pp is the mean density, Xr is the mass fraction of rock, tj, 70 j, and Hj are the half- 
life, initial abundance and decay energy per unit mass of radioisotope j, and T is the mean 
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temperature. The conductivity k{T) and specific heat Cp{T) are given in terms of their 
respective values for rock (r) and ice (i) by (e.g., Multhaup and Spohn 2007) 



k{T) = Kfr + (1 - fr)HT) 
Cp{T) = CrXr + (1 — Xr)Ci{T) 



(9) 



where is the rock volume fraction. In addition, the conductivity and density are adjusted 
by factors of 1 — 0^/^ and 1 — 0, respectively, to account for porosity (Smoluchowski 1981). 
Equation 8 yields the temperature of the planetesimal subject to radionuclide heating and 
conductive cooling. The relevant parameters for these calculations arc given in Table 2. 

A variety of outcomes are possible depending on poorly constrained model parameters, 
but melting is likely to take place if the object incorporates a significant amount of ^^Al. 
A study by Merk and Prialnik (2006) that includes radioactive heating, accretion, transfor- 
mation of amorphous to crystalline ice and melting of water ice in the formation of trans- 
neptunian objects with sizes ~ 10 km finds that the occurrence of liquid water may be 
common. Since Kuiper belt objects form in longer timescale and accrete in a colder en- 
vironment, it is not surprising that planetesimals forming early-on in the Jupiter-Saturn 
region are prone to melting. In Fig. 7 we show the mean temperature as function of time 
for several values of the delay between CAIs and planetesimal formation (assuming ^^Al is 
homogeneously distributed out to Saturn). The nebula temperature is 100 K. The values for 
the physical parameters used in this model are given in Table 1. We start with crystallized 
ice (and no CO), due to our higher temperature of accretion. We ignore serpentinization 
reactions, which would lead to further heating. Melting of these objects is difficult to avoid 
in the Jupiter-Saturn region even with significant delays in the planetesimal formation time. 
The reason is that the energy supply of ^^Al greatly exceeds the amount of energy required 
for melting the entire object. For a planetesimal of size ~ 10 km, conduction cannot remove 
this energy in time, leading to at least partial differentiation. This also fits with more de- 
tailed studies of hydrological activity of small icy planetesimals in the asteroid belt (Young 
2001), where the background temperature is higher. The disruption of partially differen- 
tiated planetesimals would then lead to a population of icy/rocky fragments available to 
ablate through the extended Kronian gas disk. 

Above we neglect heat transport by water vapor flow through a porous medium driven 
by internal subhmation (Prialnik and Podolak 1995). Yet, the presence of liquid water is not 
required for partial differentiation to take place. Enhanced conductivity by water vapor flow 
through porous media may discourage melting, but would still result in a layered internal 
state as the core would become water-depleted (Prialnik et al. 2008). To check this we we 
add a cooling term due to the rate of vapor sublimation given by (e.g., Prialnik and Podolak 



where yS"^ = — (3/ro)0 In is the surface to volume ratio (Koponen et al. 1997), ro is 
the grain size, Hy is the latent heat of vaporization, V{T) is the saturated vapor pressure. 



1995) 




(10) 
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r = 10 km, = 0.7 




time (Myr) 

Figure 7: Mean planetesimal temperature as a function of time. The curves correspond to different 
delay times between CAIs and the formation time of the primordial 10 km planetesimal. The rock 
fraction is Xr = 0.7 by mass. The porosity is (j) = 0.3. Since planetesimals in the Jupiter-Saturn 
region form in the first ~ 1 Myr they are expected to reach melting temperatures. 

Py is the gas pressure, /Xjce is the molecular mass of ice, and Rg is the gas constant. The 
value for the porosity is uncertain and the sublimation rate is model-dependent. We treat 
it as an adjustable parameter and find that for values relevant to the Jupiter-Saturn region 
with nebula temperature ~ 100 K and significant quantities of ^^Al one must sublimate a 
significant fraction of the water-ice of the object in order to keep it from melting. Thus, 
whether sublimation is sufficiently effective to preclude melting or not, radiogenic heating 
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is likely to result in separation of rock and water in our regime of interest where nebula 
temperatures are higher and accretion timescales shorter than in the Kuiper belt, which 
leads us to expect that the primordial population of planetesimals in the Jupiter-Saturn 
region is partially differentiated. 

Table 2: Thermal evolution parameters 
Thermal conductivity of ice" 

ki = 4.8812 X 107 T + 0.4685 x 10^ (g cm s'^) 

Thermal conductivity of water 

kyj = 0.55 X 10^ (g cm s"^) 

Thermal conductivity of rock** 

kr = 4.2 X 10^ (g cm s^^^) 

Specific heat of ice" 

a = 7.037 X 10^ r+ 1.85 X 10^ (cm2 K'^ s-^) 
Specific heat of water 

= 4.19 X 10^ (cm^ s-^) 
Specific heat of rock** 

= 1.2 X 10^ (cm^ s-2) 
Half life of ^^Al = 7.2 x 10^ yrs 

Initial abundance of ^^Al = 7.0 x 10"^ from ^ = 5 x 10"^ 
Heat production rate'' of ^^Al = 1.16 x 10^''' (erg g~^) 

"Hobbs (1974); ''Ellsworth and Schubert (1983); ^Castillo-Rogez et al. (2009) 



5 Conclusions 

There are two self-consistent models of satellite formation (MEa,b; Estrada and Mosqueira 
2006; see Estrada et al. 2009 for a review). Both of these models treat planetesimal dynamics 
explicitly and derive the sohds leading to the formation of the regular sateUites from the 
"shower" of planetesimal fragments following the formation of the giant planets. However, 
only one of the two accretes satellites in the presence of a relatively massive (SEMM) gas 
component (MEa,b; see Mosqueira et al. 2009 for a review) able to ablate and capture a 
significant fraction of the planetesimal rubble crossing the circumplanetary disk. Simply put, 
there are two basic options available to modelers: first, guided by the distribution of sateUite 
mass for Jupiter and Saturn (especially but not exclusively the large separation between 
Titan and lapetus), one can infer a gaseous subnebula with a relatively dense inner region 
and a long tail out to the location of the irregular satellites; second, assuming instead that 
the gas component is dissipated by an unknown mechanism in a short timescale, one can 
consider a gas-poor satellite formation model. 

Here we show that it is possible to understand the mass budgets and compositions of 
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the regular satellites of the giant planets in terms of a two-component gaseous model^. In 
particular, the result that (even allowing for porosity) lapetus is water-ice rich compared 
to solar mixtures provides strong support for the presence of a two-component (SEMM) 
gas disk at the time of satellite formation. Indeed, a lower surface density outer disk is 
required to simultaneously satisfy the constraints imposed by lapetus' small mass and ice- 
rich composition (compared to Callisto, Ganymede and Titan). Thus, lapetus' distance 
from the primary sets it apart from all other regular sateUites. But we stress that our 
framework applies in general to the formation of regular satellites around gaseous giants. In 
fact, our results also account for the water-ice enrichment of Callisto compared to Ganymede 
as observed. 

For our model to apply, the first generation of planetesimals in the Jupiter-Saturn re- 
gion must be at least partially differentiated. We conduct a preliminary investigation of the 
thermal evolution of planetesimals, and find that 10 km bodies in Jupiter-Saturn region are 
more likely to melt and differentiate than similar composition and size objects located in 
the Kuiper belt region because of the shorter accretion times and hotter nebula environment 
closer to the Sun. Since we argue that the first generation of planetesimals in the Jupiter- 
Saturn region form quickly and incorporate significant quantities of ^^Al, these objects are 
at least partially differentiated. As a result, the coUisional cascade following giant planet 
formation leads to the production of a non-homogeneous population of planetesimal frag- 
ments available to ablate through the circumplanetary disks of the giant planets, resulting in 
water-ice enrichment in the outer disk. Given the small degree of radial mixing expected in 
the quiescent- disk SEMM satellite formation model (MEa,b), a radial compositional gradient 
can be preserved in the bulk properties of the regular sateUites as observed. 

We do not include the abundances of other volatile species, as we do not have good 
constraints for their contribution to the bulk compositions of either Kuiper belt objects or 
regular satellites. Also, we do not consider a phase of condensed, refractory carbon (in 
contrast to CH4, which is trapped in water-ice). Solid carbon affects the condensate density 
and reduces the amount of CO available to sequester oxygen (Johnson and Lunine 2005; 
Wong et al. 2008). However, there is presently no evidence that such a phase plays a 
significant role in the bulk compositions of the regular sateUites, and no reason for including 
it in lapetus but not doing so in other regular satellites, Phoebe, or Kuiper belt objects. 
Although a number of the parameters employed here are at present uncertain, our work ties 
regular satellite formation to the outer solar nebula, so that constraints from observational 
studies of the compositions of the giant planets, or comet populations, and further theoretical 
work of the requirements for giant planet core formation, or Kuiper belt constraints, can be 
used in connection with the regular satellites as well. In fact, this is the first study to 
establish a direct link between solar nebula planetesimals and the bulk properties of regular 
satellites, thus opening up an entirely new research direction. A number of effects have been 
left for further work, notably: coUisional damping in the planetesimal dynamics following 
giant planet formation, and the effect of multiple disk crossings for those objects that are 

two component model has the added benefit of providing a natural explanation for the Callisto- 
Ganymede dichotomy without resorting to fine-tuning model parameters (see Estrada et al. 2009 and 
Mosqueira et al. 2009). This also fits with the observation that the composition of Callisto is icier than that 
of Ganymede. 
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not captured during the first pass. Also, our results assume that the size distribution and 
size cut-offs, and velocity at infinity of rocky and icy fragments are the same. 

However, we have tested parameter sensitivities and find that the result that the outer 
disk of Saturn is ice-rich compared to its inner disk, and also compared to the outer solar 
nebula, is robust, provided that primordial planctcsimals are at least partially differentiated. 
The reason is that delivery of solids in the outer disk is dominated by melting and vaporiza- 
tion of water-ice at temperatures of hundreds of degrees Kelvin. This result is a consequence 
of the lower gas surface density and lower planetesimal velocities in the outer disk. Conse- 
quently, icy fragments of size ^ 1 m are delivered by ablation (and some capture), whereas 
delivery of rocky fragments of similar size is highly inefficient. This is mainly because the 
temperature of objects crossing the outer disk is not high enough for rock to melt or vaporize. 
Additionally, overall a higher proportion of the rock content is delivered by capture, which 
drops off as the surface density gets lower. This is because for less pronounced ablation the 
cross section of the traversing planetesimal and gas drag force are larger, resulting in a higher 
proportion of mass captured compared to ablated. Also, rocky objects are more difficult to 
capture for a given object size due to their greater density but this may be compensated by 
the availability of somewhat smaller rocky fragments. 

Still, much work remains to be done to pin down the degree of enrichment quantitatively; 
in particular, while second pass-capture of rocky fragments does not erase a subnebula com- 
positional gradient signature, it can deliver rock to the disk via capture and should be 
explored. Also, while collisional damping of planetesimal eccentricities does not prevent 
delivery of collisional rubble to the circumplanetary disk during close encounters with the 
giant planet, the amount of mass delivered is likely to be affected, as it depends on the mass 
spectrum. Moreover, other delivery mechanisms, such as collisional capture and three-body 
effects, may also contribute to the mass budget of the regular sateUites but have not been 
included in our analysis. Finally, we assume a starting condition of a gas disk without solid 
content. Yet planetesimals dissolve in the gas giant envelope, thus enhancing the planet's 
metallicity (Pollack et al. 1986) and possibly the gas disk as well. Most of this mass delivery 
takes place before the mass of the gaseous envelope grows larger than that of the core. This 
is then followed by a dilution during runaway gas accretion (Pollack et al. 1996). It is unclear 
whether this source of solids will contribute significantly to the regular satellites, especially 
to those located far from the planet. We expect that the outer disk forms after envelope 
collapse, as the giant planets open a gap in the nebula, and that this gas is relatively free of 
solids (see Estrada et al. 2009 for further discussion). Moreover, dust setthng in the envelope 
of the giant planet may also hmit the solid content that is left behind in the subnebula far 
from the planet. Ice/rock fractionation is also possible. We leave this issue for further work. 
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